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Abstract 

We discuss a particularly symmetric model of neutrino mixings where, with good 
accuracy, the atmospheric mixing angle 623 is maximal, #13 = and the solar angle 
satisfies sin 2 612 = | (Harrison-Perkins-Scott (HPS) matrix). The discrete symme- 
try A,± is a suitable symmetry group for the realization of this type of model. We 
construct a model where the HPS matrix is exactly obtained in a first approxima- 
tion without imposing ad hoc relations among parameters. The crucial issue of the 
required VEV alignment in the scalar sector is discussed and we present a natural 
solution of this problem based on a formulation with extra dimensions. We study 
the corrections from higher dimensionality operators allowed by the symmetries of 
the model and discuss the conditions on the cut-off scales and the VEVs in order for 
these corrections to be completely under control. Finally, the observed hierarchy of 
charged lepton masses is obtained by assuming a larger flavour symmetry. We also 
show that, under general conditions, a maximal 623 can never arise from an exact 
flavour symmetry. 
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1 Introduction 



By now there is convincing evidence for solar and atmospheric neutrino oscillations. The 
Am 2 values and mixing angles are known with fair accuracy £Q. For Am 2 we have: 
Am 2 4m ~ 2.5 eV 2 and Am 2 oZ ~ 8 eV 2 . As for the mixing angles, two are large 
and one is small. The atmospheric angle #23 is large, actually compatible with maximal 
but not necessarily so: at 3a: 0.31 < sin 2 #23 < 0.72 with central value around 0.5. The 
solar angle #12 is large, sin 2 # 12 ~ 0.3, but certainly not maximal (by about 5-6 a now 
j2]). The third angle #13 is strongly limited, mainly by the CHOOZ experiment, and has 
at present a 3a upper limit given by about sin 2 #13 < 0.08. 

In spite of this experimental progress there are still many alternative routes in con- 
structing models of neutrino masses. This variety is mostly due to the considerable am- 
biguities that remain. First of all, it is essential to know whether the LSND signal 
which has not been confirmed by KARMEN [I] and is currently being double-checked by 
MiniBoone [5], will be confirmed or will be excluded. If LSND is right we probably need at 
least four light neutrinos; if not we can do with only the three known ones, as we assume 
here in the following. As neutrino oscillations only determine mass squared differences a 
crucial missing input is the absolute scale of neutrino masses (within the existing limits 
from terrestrial experiments and cosmology |01,jZj)- Even for three neutrinos the pattern 
of the neutrino mass spectrum is still undetermined: it can be approximately degenerate, 
or of the inverse hierarchy type or normally hierarchical. Given for granted that neutrinos 
are Majorana particles, their masses can still arise either from the see-saw mechanism or 
from generic dimension-five non renormalizable operators. 

At a more direct level, we do not know how small the mixing angle #13 is and how close 
to maximal is 623- One can make a distinction between "normal" and "special" models. 
For normal models 623 is not too close to maximal and 8 13 is not too small, typically a 
small power of the self-suggesting order parameter s/r, with r = Am 2 oi /Am 2 im ~ 1/35. 
Special models are those where some symmetry or dynamical feature assures in a natural 
way the near vanishing of #13 and/or of #23 — tt/4. Normal models are conceptually more 
economical and much simpler to construct. We expect that experiment will eventually 
find that #13 is not too small and that 623 is sizably not maximal. But if, on the contrary, 
either 613 very small or 823 very close to maximal will emerge from experiment, then theory 
will need to cope with this fact. Thus it is interesting to conceive and explore dynamical 
structures that could lead to special models in a natural way. 

We want to discuss here some particularly special models where both 613 and 623 — 7r/4 
exactly vanish 1 . Then the neutrino mixing matrix Ufi (f = e,fi,r, * = 1,2,3), in the basis 
of diagonal charged leptons, is given by, apart from sign convention redefinitions: 



where Ci 2 and s i2 stand for cos6 , 12 and sin 6*12, respectively. It is much simpler to write 

^^More precisely, they vanish in a suitable limit, with correction terms that can be made negligibly 
small. 
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natural models of this type with s 12 small and thus many such attempts are present in the 
early literature. More recently, given the experimental value of #12, the more complicated 
case of S12 large was also attacked, using non abelian symmetries, either continuous or 
discrete jH EE EH El CHI EH E] In many examples the invoked symmetries are particularly 
ad hoc and/or no sufficient attention is devoted to corrections from higher dimensional 
operators that can spoil the pattern arranged at tree level and to the highly non trivial 
vacuum alignment problems that arise if naturalness is required also at the level of vacuum 
expectation values (VEVs). 

An interesting special case of eq. (JTJ) is obtained for s 12 = l/v^3, i-e. the so-called 
tri-bimaximal or Harrison- Perkins-Scott mixing pattern (HPS) [13], with the entries in 
the second column all equal to l/y/3 in absolute value: 
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This matrix is a good approximation to present data 2 . It would be interesting to find a 
natural and appealing scheme that leads to this matrix with good accuracy. In fact this 
is a most special model where not only # 13 and #23 — vr/4 vanish but also 9\2 assumes a 
particular value. Clearly, in a natural realization of this model, a very constraining and 
predictive dynamics must be underlying. We think it is interesting to explore particular 
structures giving rise to this very special set of models in a natural way. In this case we 
have a maximum of "order" implying special values for all mixing angles: at the other 
extreme, anarchical models have been proposed [T3], where no structure at all is assumed 
in the lepton sector, so that, for example, #13 and 6*23 are predicted to be in no way special, 
except that there must be a smallest angle (probably near to the present bound) and a 
largest angle (expected sizably different from maximal). 

Interesting ideas on how to obtain the HPS mixing matrix have been discussed in refs 
|13j . The most attractive models are based on the discrete symmetry A4, which appears 
as particularly suitable for the purpose, and were presented in ref. (THl [TTj . In the present 
paper we start by discussing some general features of HPS models. We then present a new 
version of an A<± model, with (moderate) normal hierarchy, and discuss in detail all aspects 
of naturalness in this model, also considering effects beyond tree level and the problem of 
vacuum alignment. There are a number of substantial improvements in our version with 
respect to Ma in ref. [H]. First, the HPS matrix is exactly obtained in a first approxima- 
tion when higher dimensional operators are neglected, without imposing ad hoc relations 
among parameters (in ref. [TT]. the equality of b and c is not guaranteed by the symme- 
try). The observed hierarchy of charged lepton masses is obtained by assuming a larger 
flavour symmetry. The crucial issue of the required VEV alignment in the scalar sector is 
considered with special attention and we present a natural solution of this problem. We 

2 In the HPS scheme tan 2 #12 = 0.5, to be compared with the latest experimental determination [2|: 
tan 2 12 =O.45±° o ;°t 
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also keep the flavour scalar fields distinct from the normal Higgs bosons (a proliferation 
of Higgs doublets is disfavoured by coupling unification) and singlets under the Standard 
Model gauge group. Last not least, we study the corrections from higher dimensionality 
operators allowed by the symmetries of the model and discuss the conditions on the cut-off 
scales and the VEVs in order for these corrections to be completely under control. 

2 General Considerations 

The HPS mixing matrix implies that in a basis where charged lepton masses are diagonal 
the effective neutrino mass matrix is given by m v = {7#psdiag(mi, m 2 , m 3 )U^ PS : 



The eigenvalues of m v are m 1; m 2 , m 3 with eigenvectors (—2, 1, 1)/a/6, (1, 1, 1)/^ and 
(0, 1, — 1)/V2, respectively. In general, apart from phases, there are six parameters in 
a real symmetric matrix like m v : here only three are left after the values of the three 
mixing angles have been fixed a la HPS. For a hierarchical spectrum m 3 » m 2 >> mi, 
m\ ~ Aml tm , m\jm\ ~ /\m 2 sol / Am 2 atm and mi could be negligible. But also degenerate 
masses and inverse hierarchy can be reproduced: for example, by taking m 3 = —m 2 = mi 
we have a degenerate model, while for mi = — m 2 and m 3 = an inverse hierarchy case 
(stability under renormalization group running strongly prefers opposite signs for the first 
and the second eigenvalue which are related to solar oscillations and have the smallest mass 
squared splitting). From the general expression of the eigenvectors one immediately sees 
that this mass matrix, independent of the values of m^, leads to the HPS mixing matrix. 
It is a curiosity that the eigenvectors are the same as in the case of the Fritzsch-Xing (FX) 
matrix but with the roles of the first and the third ones interchanged (so that for HPS 
023 is maximal while sin 2 2#i 2 = 8/9, while for FX the two mixing angles keep the same 
values but are interchanged). 

If the atmospheric mixing angle is really maximal as in the HPS ansatz or close to 
maximal, it seems quite natural to interpret this as the effect of a flavour symmetry. It 
would be tempting to think of an approximate flavour symmetry such that 623 = 7r/4 
arises in the limit of exact symmetry, that is by neglecting all symmetry breaking effects. 
Here we will show that this is not the case and that, under quite general conditions, we 
can never obtain 823 = 7r/4 as a result of an exact flavour symmetry 3 . We assume that 
this symmetry is a meaningful symmetry, that is it is only broken by small effects, in 
the real world. In other words here we exclude symmetries that need breaking terms of 
order one to describe the observed fermion masses and mixing angles. Apart from that 
the symmetry can be of whatever type, global or local, continuous or discrete. Being 
interested in the limit of exact symmetry, we can neglect the sector giving rise to flavour 
symmetry breaking. We assume that the fields on which such symmetry acts are the fields 
of the standard model, plus possibly the right-handed neutrinos, so that our results will 

3 For related observations see ref. |17j . 
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also cover the see-saw case. Last, we assume canonical kinetic terms, so that the symmetry 
acts on the fields of the standard model through unitary transformations. 

Since the flavour symmetry is broken only by small effects, the mass matrices for 
charged leptons and neutrinos can be written as: 

m e = + ... , m u = ml + ... (4) 

where dots denote symmetry breaking effects and m° e has rank less or equal than one. 
Rank greater than one, as for instance when both the tau and the muon have non- vanishing 
masses in the symmetry limit, is clearly an unacceptable starting point, since the difference 
between the two non-vanishing masses can only be explained by large breaking effects, 
which we have excluded, or by a fine-tuning, which we wish to avoid. If the rank of m e 
vanishes, than all mixing angles in the charged lepton sector are undetermined in the 
symmetry limit and #23 is also completely undetermined. Therefore we can focus on the 
case when mP e has rank one. If m° e has rank one, then by a unitary transformations we can 
always go to a field basis where 

/ \ 

m° = (5) 
v m° T j 

As in the original basis, the action of the flavour symmetry on the new field basis is 
perfectly defined. If U v and U e are the unitary matrices that diagonalize m° v and m^m", 
it will be possible to adopt the parametrization 



u v = K v R 2Z {e^)p^R X z{e^)PR 12 {ei 2 ) , (6) 

where is the orthogonal matrix representing a rotation in the ij sector, P = diag(l, 1, exp iS) 
and K = diag(exp iai, exp ia 2 , exp ia 3 ). Moreover: 

U e = R 12 (9 e 12 ) (7) 

where the angle Q\ 2 is completely undetermined. The physical mixing matrix is Upmns — 
U}U U and we find: 



tan# 23 | = |cos^ 2 tan^e iQ2 + sin^ 2 ^W (<5+ai) | . (8) 
1 1 ~ " 12 cos 0£ 3 1 V ; 



Therefore, in general, the atmospheric mixing angle is always undetermined at the leading 
order. When small symmetry breaking terms are added to m° e and m°, it is possible to 
obtain 9 2 % = tt/4, provided these breaking terms have suitable orientations in the flavour 
space. If the breaking terms are produced by a spontaneous symmetry breaking through 
the minimization of the potential energy of the theory, in general two independent scalar 
sectors are needed. One of them communicates the breaking to charged fermions and the 
other one feeds the breaking to neutrinos. In such a framework a maximal atmospheric 
mixing angle is always the result of a special vacuum alignment. 

In the literature there are symmetries predicting 8 23 large, not necessarily maximal, in 
the limit of exact symmetry [T^]. For instance, this is produced by U(l) flavour symme- 
tries, when the U(l) charges of left-handed leptons and right-handed charged leptons are 
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(gi, 0, 0) and (pi,P2,0), respectively, with q 1 and pi^ all non-vanishing and different. In 
the symmetry limit, such an assignment implies (m e ~ RL, m v ~ L T L): 




( \ 

\a\ 2 a/3 
\0 a? |/3| 2 ) 



(9) 



and: 



/ \ 

m°Jm° u = |a'| 2 + |y| 2 a 7 7 ' + Y/3 / 
v a'^ + i'p \P'\ 2 + \i\ 2 j 



(10) 



with a, (3, a', j3' and 7' independent parameters of the same order of magnitude. If there 
is no conspiracy among these parameters, the resulting 823 mixing is generically large. 

In conclusion, a large lepton mixing in the 23 sector is possible as the result of an 
exact flavour symmetry. But if we want to reproduce 823 = vr/4 in some limit of our 
theory, necessarily this limit cannot correspond to an exact symmetry in flavour space. A 
maximal atmospheric mixing angle can only originate from breaking effects as a solution 
of a vacuum alignment problem. 

3 Basic Structure of the Model 

Our model is based on the discrete group A4 following refs [TU1 HI] , where its structure 
and representations are described in detail. Here we simply recall that A4 is the discrete 
symmetry group of the rotations that leave a tethraedron invariant, or the group of the even 
permutations of 4 objects. It has 12 elements and 4 inequivalent irreducible representations 
denoted 1, 1', 1" and 3 in terms of their respective dimensions. Introducing u, the cubic 
root of unity, u = exp i 2 ^-, so that 1+u+u 2 = 0, the three one-dimensional representations 
are obtained by dividing the 12 elements of A 4 in three classes, which are determined by 
the multiplication rule, and assigning to (class 1, class 2, class 3) a factor (1, 1, 1) for 1, or 
(1, u, uo 2 ) for 1' or (1, u 2 , u) for 1". The product of two 3 gives 3 x 3 = 1 + V + 1" + 3 + 3. 
Also 1' x V = 1", V x 1" = 1, 1" x 1" = V etc. For 3 ~ ( ai ,a 2 ,a 3 ), 3' ~ (61,63,63) the 
irreducible representations obtained from their product are: 



Following ref. (TTj we assigns leptons to the four inequivalent representations of A 4 . left- 
handed lepton doublets / transform as a triplet 3, while the right-handed charged leptons 
e c , /i c and r c transform as 1, 1' and 1", respectively. The flavour symmetry is broken by 
two real triplets ip and if' and by a real singlet £. At variance with the choice made by 



1 = aibi + a 2 b 2 + CI363 



(11) 



1' = a x bi + Lua 2 b 2 + uj 2 a 3 b 3 
1" = a\bx + uj 2 a 2 b2 + uja 3 b 3 

3 ~ (0263,0361,0162) 
3 ~ (0362,0163,0261) 



(12) 
(13) 
(14) 
(15) 
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jTTj, these fields are gauge singlets. Hence we only need two Higgs doublets h u>d (not three 
generations of them as in ref. fl]), which we take invariant under v4 4 . We assume that 
some mechanism produces and maintains the hierarchy (h u>d ) = v u>d <C A where A is the 
cut-off scale of the theory 4 . The Yukawa interactions in the lepton sector read: 



C Y = y e e c (pl) + y^ c {pl)" + y T r c ^iy + x a £(ll) + x d {<p'll) + h.c. + 



(16) 



In our notation, (33) transforms as 1, (33)' transforms as 1' and (33)" transforms as 1". 
Also, to keep our notation compact, we use a two-component notation for the fermion 
fields and we set to 1 the Higgs fields h Ujd and the cut-off scale A. For instance y e e c (tpl) 
stands for y e e c (pl)h d / A, x a t;(U) stands for x a ^{lh u lh u ) / A 2 and so on. The Lagrangian 
Ly contains the lowest order operators in an expansion in powers of 1/A. Dots stand for 
higher dimensional operators that will be discussed in section 6. Some terms allowed by 
the flavour symmetry, such as the terms obtained by the exchange ip' «-> ip, or the term 
ill) are missing in Ly. Their absence is crucial and will be motivated later on. 

As we will demonstrate in section 5, the fields <p', <p and £ develop a VEV along the 
directions: 



{if/) = (t/,0,0) 
(ip) = (v,v,v) 

<£> = «. 



(17) 



Therefore, at the leading order of the 1/A expansion, the mass matrices m; and m v for 
charged leptons and neutrinos are given by: 



/ 



mi = v d 



A 



Ve Ve 



Ve 



y» y^ y^ 
V Vr y^ 2 y T u J 



where 



a = x a 



K 
A 



u 
A 



/ a \ 

a d 
i d a J 

d = x d 



A 



Charged leptons are diagonalized by 



\/3 



and charged fermion masses are given by: 



(I 
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1 
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18) 



(19) 



(20) 



(21) 



VSy ( 



v d 



A 



in 



V^>y^v d 



A 



m 



T = V3y T 



Vd 



A 



(22) 



4 This is the well known hierarchy problem that can be solved, for instance, by realizing a supersym- 
metric version of this model. 
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We can easily obtain a natural hierarchy among m e , m M and m r by introducing an addi- 
tional U(l)i? flavour symmetry under which only the right-handed lepton sector is charged. 
We assign F-charges 0, 2 and 3-^4 to r c , /i c and e c , respectively. By assuming that a flavon 
9, carrying a negative unit of F, acquires a VEV (0)/A = A < 1, the Yukawa couplings 
become field dependent quantities y e ^ T = y e ,^,r{0) and we have 



y T nO{l) 



V, « 0(A 2 



ye « 0(A 



3-^4 \ 



(23) 



In the flavour basis the neutrino mass matrix reads 



m,. 



A 



/ a + 2d/3 -d/3 -d/3 \ 

-d/3 2d/3 a -d/3 
i —d/3 a — d/3 2d/3 j 



and is diagonalized by the transformation: 



U T m u U = -^diag(a + d, a, —a + d) 



with 



U 



( ^2/3 1/V3 \ 

-1/V6 l/y/3 -1/V2 
[ -1/V6 l/v^ +1/V2 ) 



(24) 



(25) 



(26) 



The leading order predictions are tan 2 #23 — 1> tan 2 ^ 12 = 0.5 and #13 = 0. The neutrino 
masses are m x = a + d, m 2 = a and m 3 = — a + d, in units of v\j 'A. We can express |a|, \d\ 
in terms of r = Am 2 oi /Am 2 tm = (|m 2 | 2 — |mi| 2 )/|m 3 | 2 — |mi| 2 ), Am 2 tm = |ms| 2 — |mi| 2 
and cos A, A being the phase difference between the complex numbers a and d: 



V2\a\ V -f 
A 

V2\d\ V f 



^ m atm 



2 cos AVI - 2r 



Vl - 2r x /Am 2 



at in 



(27) 



To satisfy these relations a moderate tuning is needed in our model. Due to the absence 
of (II) in eq. (JT5j) which we will motivate in the next section, a and d are of the same order 
in 1/A, see eq. (120)) . Therefore we expect that |a| and |d| are close to each other and, to 
satisfy eqs. ()27j) . cos A should be negative and of order one. We obtain: 



mi 



"72 



— r + 



8 cos 2 A(l - 2r) 



Am 



atm 



8 cos 2 A(l - 2r) 



Am 



atm 



1 - r + 



8 cos 2 A(l - 2r) 



Am 



atm 



(28) 



5 Notice that a unitary change of basis like the one in eq. (|21|l will in general change the relative phases 
of the eigenvalues of m„ . 
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If cos A = 



— 1, we have a neutrino spectrum close to hierarchical: 



\ma\m 0.053 eV , \rrn\ » \m 2 \ » 0.017 eV . (29) 

In this case the sum of neutrino masses is about 0.087 eV. If cos A is accidentally small, the 
neutrino spectrum becomes degenerate. The value of \m ee \, the parameter characterizing 
the violation of total lepton number in neutrinoless double beta decay, is given by: 



\m e 



1 + 4r ^ 



9 8cos 2 A(l - 2r) 



Ami • (30) 



For cos A = — 1 we get \m ee \ ~ 0.005 eV, at the upper edge of the range allowed for normal 
hierarchy, but unfortunately too small to be detected in a near future. Independently from 
the value of the unknown phase A we get the relation: 

2 2 10 2 / r\ 

\m 3 \ =\m ee \ + —km atm \l- -1 , (31) 

which is a prediction of our model. 

It is also important to get some constraint on the mass scales involved in our construc- 
tion. From eqs. (|27jl and ()20j) . by assuming Xd ~ 1 v u ~ 250 GeV, we have 

A « 1.8 X 10 15 f y GeV . (32) 

Since, to have a meaningful expansion, we expect v ' < A, we have the upper bound 

A < 1.8 x 10 15 GeV . (33) 

Beyond this energy scale, new physics should come into play. The smaller the ratio v'/A, 
the smaller becomes the cut-off scale. For instance, when v'/A = 0.03, A should be close to 
10 14 GeV. A complementary information comes from the charged lepton sector, eq. (|22j) . 
A lower bound on v/A can be derived from the requirement that the Yukawa coupling y T 
remains in a perturbative regime. By asking y T Vd < 250 GeV, we get 

v > 0.004 . (34) 
A 

Finally, by assuming that all the VEVs fall in approximately the same range, which will 
be shown in section 5, we obtain the range 

v v 11 

0.004 < — « - « - < 1 , (35) 
AAA v ' 

that will be useful to estimate the effects of higher-dimensional operators in section 6. 
Correspondingly the cut-off scale will range between about 10 13 and 1.8 x 10 15 GeV. 



8 



4 Vacuum alignment 



In this section we investigate the problem of achieving the vacuum alignment of eq. (|17|) . 
At the same time we should prevent, at least at some level, the interchange between 
the fields ip and ip' to produce the desired mass matrices in the neutrino and charged 
lepton sectors. As we will see, there are several difficulties to naturally accomplish these 
requirements. By minimizing the scalar potential of the theory with respect to (p and ip' 
we get six equations that we would like to satisfy in terms of the two unknown v and v'. 
Even though we expect that, due to the symmetry A 4 , the six minimum conditions are 
not necessarily independent, such an expectation turns out to be wrong in the specific 
case, unless some additional relation is enforced on the parameters of the scalar potential. 
These additional relations are in general not natural. For instance, even by imposing 
them at the tree level, they are expected to be violated at the one-loop order. Therefore, 
as we will now illustrate, the minimum conditions cannot be all satisfied by our vacuum 
configuration. 

As an example here we analyze the most general renormalizable scalar potential invari- 
ant under A4 and depending upon the triplets ip and p' of the Lagrangian C in eq. 
The term (11) in L can be forbidden by an additional symmetry, commuting with A*. One 
possibility is just the total lepton number L or a discrete subgroup of it. Here we consider 
a Z4 symmetry under which f c transform into — if c (/ = e,/i, r), / into il, ip is invariant 
and ip' changes sign. This symmetry also explains why <p and <p' cannot be interchanged. 
The scalar potential V contains bilinears Bi, trilinears T, and quartic terms Qi, invariant 
under the group At x Z 4 . A choice of independent invariants is: 



B 2 



2 1 2 1 2 

p x + p 2 + <^3 

' 2 1 / 2 / 
<Pl +^2 +^3 



1 



ip X p' 2 p' z + P2P Z P' X + <p3<pWs 



Qi 
Q2 

Q 4 
Q 5 

Q 7 



<pI<p 2 2 + pipi + p 2 3 pi 
1 2 1 2 2 1 2 1 2 

iPi+uJ ip 2 + up 3 \ 

I 2 1 2 / 2 / 2 1 2 1 2 

tp x p 2 + p 2 ip 3 + ip 3 p 1 

I / 2 . 2 /2 . / 2 ,2 

\<Pl ip 2 +UP 3 I 

p 1 p 2 p\p' 2 + p 2 p 3 p' 2 <p' 3 + (p3<Pl<fh<f{ 

(p\ + pj + <pi){p[ 2 + p'2 + y>' 3 2 ) 

(pi + u?p\ + ujpl)(p[ 2 + ujp' 2 2 + uVf) 



(36) 



The scalar potential reads: 



V 



+ c\Q\ + C2Q2 + C3Q3 + C4Q4 

+ C5Q5 + caQa + (c 7 Q 7 + c.c) 



(37) 
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We start by analyzing the field configuration: 



((p) = (v,v,v) 

(<p') = (v',0,0) (38) 

The minimum conditions are: 

M 2 v + ii X v 2 + 4dv 3 + 2c 6 vv' 2 + 2(c 7 + c 7 )vv' 2 = 
M\v + ^iv 2 + 4dv 3 + 2c 6 vv' 2 + 2(uj 2 c 7 + uc 7 )vv' 2 = 
M\v + /iiw 2 + Aciv 3 + 2c 6 vv' 2 + 2(ujc 7 + u 2 c 7 )vv' 2 = 
M 2 V + 4c 4 v' 3 + 6c 6 v 2 v' = 
H2vv' + c^v 2 v' = 

fi 2 w' + c 5 v 2 v' = . (39) 

The equations dV/d(fi = are clearly incompatible unless c 7 = 0. Even by forcing c 7 to 
vanish, we are left with three independent equations for the two unknown v and v', which, 
for generic values of the coefficients, admit only the trivial solution v — v' — 0. This 
negative results cannot be modified by adding to V the terms depending on the singlet £. 
Also by investigating the problem in a slightly more general framework, with (p real and 
(<//, £) complex, we reach the same conclusion. Although we have not a no-go theorem, 
these examples show the difficulty to obtain the desired alignment. 

The difficulty illustrated above is not common to all vacua. For instance the other 
possible alignment: 

(tfi) = (v,v,v) 

(<p') = (v',v',v') (40) 

leads to the minimum conditions: 

M 2 v + ^iV 2 + /2 2 v' 2 + Aciv 3 + (2c 5 + Qc 6 )vv' 2 = 

M 2 V + 2fi 2 vv' + 4c 3 v' 3 + (2c 5 + 6c 6 )v 2 v' = . (41) 

In a non- vanishing portion of the parameter space, these equations have non-trivial solution 
with non- vanishing v and v'. 

It is possible to show that, by sufficiently restricting the form of the most general scalar 
potential invariant under At, the desired alignment can be obtained. Restrictions that are 
unnatural in a generic model becomes technically natural in a supersymmetric (SUSY) 
model. The well-known non-renormalization properties of the superpotential allow to 
accept, at least from a technical viewpoint, a restricted number of terms, compared to 



dV 

d(fi 
dV 

dip 2 
dV 

d(p 3 
dV 

Wi 

dV 
dip' 2 
dV 
d<p' 3 



dV 

difi 
dV 
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what the A 4 symmetry would permit. Undesired terms of the superpotential that are set 
to zero at the tree level are not generated at any order in perturbation theory. Indeed 
we have produced a SUSY example of this type, where the alignment problem is solved 
and this example is discussed in detail in the Appendix. However our real aim is to build 
a fully natural model, where all the terms allowed by the symmetries are present and 
where the only deviations from the symmetry limit are provided by higher- dimensional 
operators, rather than by small violations of ad-hoc imposed relations. As we will now 
see, there exist a simple and economic solution in the context of theories with one extra 
spatial dimension. 

5 A4 model in an extra dimension 

One of the problems we should overcome in the search for the correct alignment is that of 
keeping neutrino and charged lepton sectors separate, including the respective symmetry 
breaking sectors. Here we show that such a separation can be achieved by means of an 
extra spatial dimension. The space-time is assumed to be five-dimensional, the product of 
the four-dimensional Minkowski space-time times an interval going from y = to y = L. 
At y = and y = L the space-time has two four- dimensional boundaries, which we will 
call branes. Our idea is that matter SU(2) singlets such as e c , /i c , r c are localized at y = 0, 
while SU(2) doublets, such as I are localized at y = L (see Fig.l). Neutrino masses arise 
from local operators at y — L. Charged lepton masses are produced by non-local effects 
involving both branes. Later on we will see how such non-local effects can arise in this 
theory. The simplest possibility is to introduce a bulk fermion, depending on all space-time 
coordinates, that interacts with e c , /i c , r c at y = and with I at y = L. The exchange of 
such a fermion can provide the desired non-local coupling between right-handed and left- 
handed ordinary fermions. Finally, assuming that <p and (ip', £) are localized respectively 
at y = and y = L, we obtain a natural separation between the two sectors. 

5.1 Alignment in an extra dimension 

Such a separation also greatly simplify the vacuum alignment problem. We can determine 
the minima of two scalar potentials Vq and Vl, depending only, respectively, on ip and 
(<£>', £). Indeed, as we shall see, there are whole regions of the parameter space where 
Vo(<p) and Vl(</?',£) have the minima given in eq. (JT7J). Notice that in the present setup 
dealing with a discrete symmetry such as A 4 provides a great advantage as far as the 
alignment problem is concerned. A continuous flavour symmetry such as, for instance, 
SO(3) would need some extra structure to achieve the desired alignment. Indeed the 
potential energy / d 4 x[Vo(<p) +Vl{<p' , £)] would be invariant under a much bigger symmetry, 
SO(3) x SO(3)i, with the SO(3) acting on ip and leaving (<p',£) invariant and vice- 
versa for SO(3)/,. This symmetry would remove any alignment between the VEVs of ip 
and those of (<//,£). If, for instance, (|17|) is minimum of the potential energy, then any 
other configuration obtained by acting on (|T7|) with SO(3)oX SO(3)l would also be a 
minimum and the relative orientation between the two sets of VEVs would be completely 
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Figure 1: Fifth dimension and localization of scalar and fermion fields. The symmetry 
breaking sector includes the triplets ip and <p', localized at the opposite ends of the 
interval. Their VEVs are dynamically aligned along the directions shown at the top of the 
figure. 



undetermined. A discrete symmetry such as A4 has not this problem, as we will show 
now. 

Consider first the scalar potential Vq (<p): 



M 2 

V {<p) = ~yB\ + ii x T x + c x Q x + c 2 Q 2 



(42) 



where Bi, T\, Qi 2 are defined in eq. (pjj) . The minimum conditions at p = (v,v,v) are: 
dV 



while the minimum condition at ip 

dV 



v{Ml + inv + Adv 2 ) = 
(v,0,0) is: 
v(M 2 + 4c 2 v 2 ) 



1,2,3) 



dip 1 







(43) 



(44) 



since in this case (dV$/ dip2$) = are automatically satisfied. Both tp = (v,v,v) and 
(p = (v, 0, 0) can be local minima of Vo, depending on the parameters. The constants 
Ci j2 should be positive, to have V bounded from below. We can look at the region where 
<^ \Mi\. When c\ 3> c 2 and M\ < 0, the minimum at ip = (v, 0, 0) is the absolute one, 
while for c 2 3> c\ and Mf < Vo is minimized by ip = (v, v, v). Therefore we have a large 
portion of the parameter space where the minimum is of the desired form: ip = (v, v, v). To 
be precise, in this region, there are four degenerate minima: ip = (v, v, v), ip = (v, —v, —v) 
<p = (—v,v, —v) p = (—v, —v,v), related by A 4 transformations. 

Now we turn to Vl{<p',£,)- As we did in section 4.1, we assume both <p' and £ real 
and odd under the action of a discrete Z4 symmetry. The most general renormalizable 
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invariant potential is a combination of B 2 , in eq. (|36|) and the following invariants: 

b 3 = e 
q 8 = e 

Q 9 = ^iVaVa 

Qio = £Vi 2 + ^> 2 + ^ 2 ) . (45) 

We have: 

M 2 M 2 

Vl (<p', = -yB 2 + -^-B 3 + C3Q3 + C4Q4 + c 8 Q s + c 9 Q 9 + c 10 Qi • (46) 
We search for a minimum at <// = (v f , 0, 0) and £ = it: 

t/(M 2 2 + 4c 4 t/ 2 + 2ci m 2 ) = 



u(M 2 + Ac 8 u 2 + 2c 10 v' 2 ) = , (47) 



while (dVi/dip^ 3 ) = are always satisfied. There is a region of the parameter space where 
the absolute minimum is of this type. Taking into account the A4 symmetry, in all this 
region we have six degenerate minima: if' = (±v 0, 0),ip' = (0, ±i/, 0) and </?' = (0, 0, ±v '). 
Putting together the minima of Vo(<p) and Vl(ip',£) we have 24 degenerate minima of the 
potential energy, differing for signs or ordering. It can be shown that these 24 minima 
produce exactly the same mass pattern discussed in section 3.1, up to field and parameter 
redefinitions. Therefore it is not restrictive to choose one of them, for instance ip = (v, v, v) 
and cp' = (v', 0, 0), to analyze the property of this model. 

The observed hierarchy among lepton masses can be efficiently described by an ad- 
ditional U(1)f flavour symmetry, under which only right-handed charged leptons are 
charged: F(e c ,/i c , r c ) = (4,2,0). To spontaneously break this symmetry and to produce 
the desired hierarchy, we need a scalar field 6, carrying a negative unit of F and developing 
a VEV {0)/A m 0.22. In our framework 9 is localized on the brane at y = and the scalar 
potential V of eq. is modified into: 

V - V + M%B 4 + Cll Qu + c 12 Q 12 , (48) 

where 

B, = \9\ 2 
Qu = \9\ 4 

Q 12 = \6\ 2 { v { + vl + v l) ■ (49) 

The minimum conditions at ip — (v,v,v) and |^| = t read: 
dV 

— ° = V (M 2 + fj, lV + A Cl v 2 + 2c 12 t 2 ) = (i = l,2,3) , 

o<pi 

dv 

= 2t(M 2 + 2 Cll t 2 + 3c 12 v 2 ) = . (50) 

These conditions are satisfied by non- vanishing (t, v) in a finite portion of the parameter 
space. Therefore the inclusion of an abelian flavour symmetry is fully compatible with the 
mechanism for vacuum alignment discussed above. 
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5.2 Lepton masses and mixing angles 



We now show how it is possible to take advantage of above results to obtain the desired 
lepton masses. To this purpose we introduce a bulk fermion field F(x, y) = (F\, F 2 ), singlet 
under SU(2) with hypercharge Y — — 1 and transforming as a triplet of A*. We also impose 
the discrete Z 4 symmetry introduced in section 4.1 under which (f c , I, F, <p, <p',£) transform 
into (—if c , il, i>F, <p, —</?', — £)• The action is 



S 



d 4 xdy 



iF 1 a»d^F l + iF 2 a^F 2 + \{F 2 d y F x - d y F 2 F x + h.c. 



- M{F X F 2 + F X F 2 ) 

+ V (<p)6(y) + V L (<p',05(y-L) 

+ [Y e e c {ipFi) + YpfjffrFj' + Y^Fj)' + h.c] 5(y) 



+ 



^{ll)h u h u + ^'U)h u h u + Y L (F 2 l)h d + /i.e.] S(y — L)| + ... 



(51) 



where the constants Y have mass dimension -1/2. The first two lines represent the five- 
dimensional kinetic and mass terms of the bulk field F. The third line is the scalar 
potential and the remaining terms are the lowest order invariant operators localized at the 
two branes. Dots stand for the kinetic terms of f c , I, ip, ip', £ and for higher- dimensional 
operators, which will be classified in section 6. 

The potential energy is given, at lowest order by: 



U 



d 4 x [V (<p) + V L (<p',0] 



(52) 



and, under the conditions discussed above, is minimized by eqs. (fTTj) . It is clear that tp 
and (<£>',£) are strictly separated only at lowest order. Indeed higher-dimensional brane 
interactions like for instance (<pipFiF 2 ) / 'A 2 , ((p'ip'FiF 2 )/A 2 are allowed. At the one-loop 
level, the exchange of the bulk fermion F will give rise to the structures (55,6,7 of eq. 
(J36)) and this will necessarily deform the vacuum (|17|) . Here we will assume that such a 
deformation is sufficiently small. Indeed, as we shall see in section 6, the operators of the 
type (ipipip'ip 1 ) arising from one-loop F-exchange, are suppressed by 1/A 4 L 4 . 

We now discuss the effects of the tree-level exchange of F. To this purpose we consider 
the equations of motion for (F 1; F 2 ): 



ia»d II F 2 + d y F 1 -MF 1 = 
ia^duFi — d v F 2 — MF 2 = 



(53) 



If M is large and positive, we can prove that all the modes contained in (Fi,F 2 ) become 
heavy, at a scale greater than or comparable to l/L, which we assume to be much higher 
than the electroweak scale. If we are only interested in energies much lower than l/L, 
we can solve the equations of motion in the static approximation, by neglecting the four- 
dimensional kinetic term: 



Fx{y) 
F 2 (y) 



F x {L)e M ^ 
F 2 (0)e- My . 



(54) 
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These equations must be supplemented with appropriate boundary conditions, which we 
can identify by varying the action S with respect the fields (Fx, F 2 ). The boundary terms 
read 



(8S) boundary = J d*X UFx(L) 



+ SF 2 (L) 
+ 6Fx(0) 
+ 5F 2 (0) 



■-Fx(L) + Y L lh d 



--F 2 (0) + Y e e e <p + + Y T r c <p T 



(55) 



where (p^ = (</>!, cjy? 2 , i^ 2 ^) and tp T = (<^i, cu 2 ^, ^3)- We can chose as boundary condi- 
tions: 



Fx(L) = 2Y L lh d 

F 2 (0) = 2(Y e e c <p + Y^ c ( p tM + Y T r c <f T ) 



(56) 



Since 5Fx(L) = <5F 2 (0) = 0, we have (tfS) boundary = 0, as desired. By substituting back 
eqs. (J54*Jl and ()56|) into the action S we get 



S = U + 



d x 



y. 



Vt 



+ ... 



with 



Vj_ 
A 



4Y L Y f e 



-ML 



(/ = e,n,r) 



(57) 



(5? 



Therefore, in lowest order approximation we have reproduced the Lagrangian Ly of eq. 
(jTUj) and the discussion of section 3 applies. 

We also recall that, to account for the observed hierarchy of the charged lepton masses, 
we have included an additional U(l) flavour symmetry. Therefore, in the present picture, 
the quantities Y e ^ r stand for: 

2 



y. = Y.(l 



Y T = Y T 



(59) 



where Y e ^ T are field-independent constants having similar values. After spontaneous 
breaking of U(l), the Yukawa couplings yj possess the desired hierarchy. 



6 Higher-order corrections 

The results of the previous section hold to first approximation. Higher-dimensional op- 
erators, suppressed by additional powers of the cut-off A, can be added to the leading 
terms in eqs. ()42I46I52I57I58|) . Here we will classify these terms and analyze their physical 
effects. In particular we will show that these corrections are completely under control in 
our model and that they can be made negligibly small without any fine-tuning. We can 
order higher-order operators into three groups. 
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6.1 Local corrections to m, 



There are higher-order operators that are local in the five-dimensional theory and do not 
depend upon the heavy fermion sector (Fi,Fz). As we have seen, at leading order, the 
neutrino mass matrix m v arises entirely from operators of this type that are localized at 
y = L. On this brane we only have scalar fields (£, (p 1 ), odd under Z4. Therefore higher- 
dimensional operators modifying m u and localized at y = L are down by two powers of 
the cut-off, compared to the leading ones. After A4 breaking, the only two operators that 
cannot be absorbed by a redefinition of the parameters x a ,d are: 



A 4 



x 



i{ip'ip)'{ll)"h u h v 



(60) 



^ v ' v ')'\U)'h u h u 



After adding these operators localized at y = L to the five- dimensional action of eq. (foT|). 
we get a neutrino mass matrix 



A 



( a+b+c \ 

a + ub + u 2 c d 
d a + uj 2 b + uc j 



where 



b = x b 



uv 
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uv 



A 3 



C — X r 



A 3 



(61) 



(62) 



to be compared with a and d of eq. (J^UJ) . 



6.2 Corrections from tree-level F-exchange 

Another set of higher-dimensional operators arise from the exchange of the heavy fermion 
(Fi, F2) in the static limit and in the tree-level approximation. To classify them, we should 
list all operators localized at the two branes that are linear in the bulk fermion (^,^2). 
At y = such operators have the generic structure 

Y f r<pFi , -j^fVFi , -j^fVFi ... (/ = e,/i,r) (63) 

After spontaneous A4 breaking, the effect of these operators can be absorbed by redefin- 
ing the coupling constants Y), (/ = e, /x, r), at least up to order ip 3 . Thus the leading 
interactions between f c and F\ 

[r B e c (pFi) + YpH c {ipFi)" + Y T T c {<pF x )' + h.c] 5(y) (64) 

are unchanged up to relative order 1/A 2 . We are left with the couplings of F 2 at the brane 
y = L. Neglecting all operators that, after v4 4 breaking, only lead to a renormalization of 
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the parameter Yl, we find four new terms: 



^W)'(F 2 l)"h d 



jit [<p'i(F 2 ) 2 l 3 + <f/ 2 {F 2 ) 3 h + ^{F 2 ) x l 2 \ h a 



(65) 



^ Wi{F 2 ) 3 l 2 + V ' 2 {F 2 ) x l 3 + ^{F 2 ) 2 h\ h d 

After the breaking of A±, the leading order interaction of F 2 at y = L is modified by the 
operators (J65|) to 

Y L (F 2 l)h d + h.c] 8{y-L) 



(66) 



where 







( l + z x + z 2 





\ 


(h 


\ 




h 







1 + uzi + oj 2 z 2 




h 




(67) 


V & ) 




V 


Z3 


1 + CJ 2 2i + UJZ 2 1 




J 





#1,2 



Zt, 2 v' 2 



z 3,4 



Z 3 4 UV ' 



(6* 



After integrating out the heavy modes in (F 1 ,F 2 ) in the limit of vanishing external mo- 
menta for the light modes, we obtain the effective four-dimensional Lagrangian 



V j-e c {ipl)h d + + ^r c ^l)'h d 



A 



AY L Y f e 



-ML 



(/ = e,/i,r) 



(69) 
(70) 



The mass matrix for the charged leptons becomes 



mi = v d — 



/ y e (l + z 1 + z 2 ) y e (l + uz 1 + u 2 z 2 + z 3 ) y e (l + u 2 z 1 + uz 2 + z 4 ) ^ 

+ z\ + z 2 ) y^uo(l + uozi+uj 2 z 2 + ujz 3 ) y^0J 2 (l + uj 2 Zi + uoz 2 + uo 2 z 4 ) 
y y T (l + z 1 + z 2 ) y T uj 2 (l + ujzi+uj 2 z 2 + uj 2 z 3 ) y T u{l + uj 2 z 1 + ujz 2 + uz 4 ) j 

(71) 



6.3 Effects on masses and mixing angles 

To first order in the small parameters b, c and Zi, the neutrino masses are modified into: 

mi = (a + d-i(6 + c))^ 
v 2 

m 2 = {a + b + c)-^ 

m 3 = {- a + d+ l -{b + c)Yj- (72) 
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and the charged lepton masses are changed into 



m P 



m. 



V3y e (l + - + T )v dJ 



V3y T (l + uj 2 j + ^yK^ 



(73) 



To the same order, but neglecting terms like Z{ y e ,^/y T , we get: 



\U, 



e3 1 



(b-c)(d-a) + (b-c)(d + a) 1 
2^/2(ad + ad) V2 



-z 1 + z 2 + -j=z 3 



Z4 



tan 2 9 2 3 1 
tan 2 6*i2 1 



(&-c)d+ (b-c)d 



(ad + ad) 



+ 2 



^2 + Z 2 + ~{Z 3 + Z 3 + Z 4 + Zi) 



«1 



_ . ^3 + ^3 . ^4 + ^4 
-22 - -22 H o 1 ~ 



(74) 



These relations explicitly show that the corrections induced by the higher- dimensional 
operators are of order uv'/A 2 or v' 2 /A 2 . From our estimate in eq. (|35|1 we see that these 
parameters can be as small as 2 10~ 5 . If the cut-off A is one order of magnitude larger that 
the VEVs of the model, the resulting corrections are at the level of one percent, already 
beyond any planned experimental test. If on the contrary, the VEVs are anomalously 
close to the cut-off A, then eqs. (|73jl show that deviations roughly of the same size are 
expected in U e3 , tan 2 #23 and tan 2 6*12. How much close to A can the VEVs be? We expect 
that the subleading corrections do not spoil the leading order form of the neutrino mass 
spectrum, eq. (|2*Hjl . This implies that v' 2 /A 2 <^ r, so that r sets the natural upper bound 
to the expected deviations from the leading order results. 



6.4 Corrections from one-loop F-exchange 

Further corrections to lepton mass matrices and to the scalar potential can arise from one- 
loop exchange of (Fx, F 2 ) in the static limit. Consider for instance the following operators 
localized at y — and at y = L: 



-<pFxF 2 5(y) , 



1 



illF l F 2 h u h u 5(y - L) 



By integrating out, at one-loop order, the heavy modes contained in (F\,F 2 ) we get: 



(75) 



■jf&UKKj d 4 kA F (k,0,L)A F (k,L,0) 



(76) 



where k is the four-momentum running in the loop and Ap(k,y,y') is the adimensional 
propagator of (Fx, F 2 ) in a mixed momentum-space representation. Since the loop integral 
is convergent, we get 

1 f(ML) 



-£ipllh u h u 



(77) 



A 3 A 4 L 4 

where f(ML) is a function of the adimensional combination ML. Thus the resulting local 
operator is suppressed by four additional powers of the cut-off scale. This behavior is quite 
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generic and similar suppressions are found for other operators originating from one-loop 
exchange of (F 1; F 2 ). 

The corrections that modify the scalar potential discussed in the previous section are 
of this type. As an example, consider the localized interactions: 

^ipipF x F 2 5{y) , ^<p'<p'FiF 2 8(y-L) . (78) 

Also in this case, after integrating over F%) in the limit of vanishing external momenta, 
we get: 

4§W • (79) 

Due to their large suppression, these corrections are negligible compared to those discussed 
above. 



7 Conclusion 

There are by now several theoretical mechanisms that can qualitatively explain the ob- 
served large lepton mixing angles [T9] . They are sufficiently flexible to quantitatively 
accommodate the measured parameters. They are also compatible with our ideas on 
quarks masses and mixing angles so that they can be nicely embedded into a unified pic- 
ture of fermion properties, such as, for instance, a grand unified theory. Many of these 
mechanisms predict a generically large atmospheric mixing angle and a generically small 
#13 angle, without favouring any specific value for these parameters. The best values of 
global fits are currently very close to 9 2 3 = tt/4 and 9\3 = 0, but the experimental errors 
still allow for large deviations from these remarkable values. Indeed, according to many 
of the above mentioned mechanisms, deviations from #23 = vr/4 and #13 = are expected 
at the observable level. It may take a long time before such deviations can be actually 
observed. A sensitivity on #13 around 0.05 is foreseen in about ten years from now, with 
the full exploitation of high-intensity neutrino beams. A reduction by a factor of two of 
the present error on 623 will also require special neutrino beams and a similar time scale. 

It might happen that after all this experimental effort, (#23 ~^/^) an d $13 still remain 
close to zero, within errors. At this point it would be legitimate to suspect that such 
special values are produced by a highly symmetric flavour dynamics. Given the already 
good experimental precision on #12, the so-called Harrison-Perkins-Scott mixing scheme, 
where 623 = vr/4, 613 = and sin 2 9 12 = 1/3, would fit very well the data. In this paper we 
have proposed a model that reproduces accurately the HPS mixing pattern. We started 
by discussing whether such a pattern can be obtained from an exact flavour symmetry. 
We showed that, under general conditions, an exactly maximal atmospheric mixing angle 
cannot arise from an exact flavour symmetry. The flavour symmetry should be necessarily 
broken and a maximal #23 is the result of a special alignment between the breaking effects 
in the neutrino sector and those occurring in the charged lepton sector. If the flavour 
symmetry is spontaneously broken, this corresponds to a non-trivial vacuum alignment. 
Our model gives rise to the HPS mixing scheme in the context of a spontaneously broken 
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A 4 flavour symmetry, A A being the discrete subgroup of SO (3) leaving a tetrahedron 
invariant. 

At leading order, that is by neglecting symmetric operators of higher dimension, neu- 
trino masses only depend on two complex Yukawa coupling constants. Due to the un- 
known phase difference between these two constants, we cannot determine the absolute 
scale of neutrino masses. We expect that the neutrino spectrum is of the normal hi- 
erarchy type but not too far from degenerate. At leading order the model predicts 
|m 3 | 2 = \m ee \ 2 + 10/9Am^ m (l — r/2). A remarkable feature of our model is that at 
the leading order the lepton mixing angles are completely independent from these two 
parameters, so that the HPS mixing pattern is always obtained. The lepton mixing de- 
pends entirely on the relative alignment between the VEVs giving masses to the neutrino 
sector and those giving masses to the lepton sector. We discuss in detail the problem of 
vacuum alignment. To avoid the proliferation of Higgs doublets, the scalar fields breaking 

are gauge singlets in our model. We propose an unconventional solution to the vac- 
uum alignment problem, where an extra dimension described by a spatial interval plays 
an important role. Two scalar sectors live at the opposite ends of the interval and their 
respective scalar potentials are minimized by the desired field configurations, for natural 
values of the implied parameters. Such a mechanism only works in the case of discrete 
symmetries, since in the continuous case the large symmetry of the total potential energy 
would make the relative orientations of the two scalar sectors undetermined. We have 
also extensively discussed how this lowest order picture is modified by the introduction of 
higher dimensional operators. The induced corrections are parametrically small, of second 
order in the expansion parameter VEV/A, A being the cut-off of the theory, and they can 
be made numerically negligible. Last but not least, the hierarchy of the charged lepton 
masses can be reproduced by the usual Froggatt-Nielsen mechanism within the context 
of an abelian flavour symmetry, which turns out to be fully compatible with the present 
scheme. 

We believe that, from a purely technical point of view, we have fulfilled our goal to 
realize a completely natural construction of the HPS mixing scheme. But to construct 
our model we had to introduce a number of special dynamical tricks (like a peculiar set of 
discrete symmetries in extra dimensions). Apparently this is the price to pay for a "special" 
model where all mixing angles are fixed to particular values. Perhaps this exercise can be 
taken as a hint that it is more plausible to expect that, in the end, experiment will select 
a "normal" model with 13 not too small and #23 not too close to maximal. 
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Appendix 



Here we discuss a SUSY solution to the vacuum alignment problem. In a supersymmetric 
context, the right-hand side of eq. (JTfij) should be interpreted as the superpotential wi of 
the theory, in the lepton sector. A key observation is that this superpotential is invariant 
not only with respect to the gauge symmetry SU(2)x U(l) and the flavour symmetry 
U(1)f x A4, but also under a discrete Z 3 symmetry and a continuous U(1)jj symmetry 
under which the fields transform as shown in table 1. 



Field 
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e c 


If 


T C 


h u ,d 


<P 
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¥0 


¥0 






3 
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1 
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z 3 


us 


u 2 


u 2 


u 2 
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us 


UJ 
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us 


UJ 


U(1)r 
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1 


1 














2 


2 


2 



We see that the Z% symmetry explains the absence of the term (11) in wf. such a term 
transforms as uo 2 under Z 3 and need to be compensated by the field £ in our construction. 
At the same time Z 3 does not allow the interchange between (p' and (p, which transform 
differently under Z3. Charged leptons and neutrinos acquire masses from two independent 
sets of fields. If the two sets of fields develop VEVs according to the alignment described 
in eq. (|17|) . then the desired mass matrices follow. 

Finally, there is a continuous U(1)r symmetry that contains the usual .R-parity as a 
subgroup. Suitably extended to the quark sector, this symmetry forbids the unwanted di- 
mension two and three terms in the superpotential that violate baryon and lepton number 
at the renormalizable level. The U(1)r symmetry allows us to classify fields into three 
sectors. There are "matter fields" such as the leptons I, e c , /i c and r c , which occur in 
the superpotential through bilinear combinations. There is a "symmetry breaking sector" 
including the higgs doublets h u d and the flavons (p 1 , <p and £. As we will see these fields 
acquire non-vanishing vacuum expectation values (VEVs) and break the symmetries of 
the model. Finally, there are "driving fields" such as <p' Q , ipo and £0 that allows to build 
a non-trivial scalar potential in the symmetry breaking sector. Since driving fields have 
R-charge equal to two, the superpotential is linear in these fields. 

The full superpotential of the model is 

w = wi + w d (80) 

where, at leading order in a 1/A expansion, wi is given by the right-hand side of eq. ()16j) 
and the "driving" term Wd reads: 

w d = M(ip (p) + g(ip <p<p) + gi(<Po<p'tp') + 92^(^>W) + SfcfofaV) + ^o^ 2 • (81) 

We notice that at the leading order there are no terms involving the Higgs fields h u ^. We 
assume that the electroweak symmetry is broken by some mechanism, such as radiative 
effects when supersymmetry (SUSY) is broken. It is interesting that at the leading order 
the electroweak scale does not mix with the potentially large scales u, v and v ' . The scalar 
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potential is given by: 



dw 



+ m^| 2 + ... (82) 



where <pi denote collectively all the scalar fields of the theory, mf are soft masses and dots 
stand for D-terms for the fields charged under the gauge group and possible additional soft 
breaking terms. Since mj are expected to be much smaller than the mass scales involved 
in Wd, it makes sense to minimize V in the supersymmetric limit and to account for soft 
breaking effects subsequently From the driving sector we have: 



dw 

dipoi 
dw 

dw 

<9y?o3 
dw 

dtp'oi 
dw 

dw 

Ms 
dw 

Wo 



M(pi + g(p 2 ^3 = 
Mip 2 + g^ipi = 
M(f 3 + gfi(f2 = 
9i<p2<Pa + 92&i = 
gicp' 3 cp[ + g 2 &2 = 
gi(p[(p2 + g2&' 3 = 

g 3 (<p'<p') + g 4 e = (83) 



The first three equations are solved by (up to irrelevant sign ambiguities): 

M . . 

(p=(v,v,v) , v — . (84) 

9 

The remaining equations are solved, in general, by: 

^' = (0,0,0) , ^ = , (85) 

unless some further relation is imposed on the coefficients gx, ...,#4. If g% = 0, then, up to 
an irrelevant reordering, we have 

if/ = (u',0,0) , 

t = u = -^W) (86) 

94 

with v' and u undetermined. In this case we find that, for m^, o ,m^/ , m| o > 0, the driving 
fields ifo, ip' Q and £ vanish at the minimum. Moreover, if mi,, mj < 0, then u and v' slide to 
large scales, eventually stabilized by one-loop radiative corrections. The supersymmetric 
case is better than the non-supersymmetric case in two respects. First of all, at least from 
a technical viewpoint, the absence of a term in the superpotential is radiatively stable. 
Moreover, as we have seen, once g 2 has been set to zero, the equations selecting (fTTj) as 
the correct minimum are consistent. 
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